Mode Coupling in Space and Time Varying 
Anisotropic Absorbing Plasmas 



Bijan Sheikholeslami Sabzevari 



Institute for Studies in Theoretical Physics and Mathematics 
P.O. Box 19395-5531, Tehran, Iran 



Abstract 



A four dimensional systematic mathematical approach for investigating propagation 
and coupling of wave modes in a slowly varying (in all space directions and time) 
anisotropic, absorbing plasma is represented. The formalism is especially useful 
for energy considerations of the waves. It is applicable to general cases of mode 
conversion in plasmas with general geometries of space-time and magnetic field con- 
figurations. A simple example of how this formalism can be applied to practical cases 
is given. 
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1 Introduction 



Almost all investigations about wave coupling and mode conversion in plasmas were done for 
cases, where the inhomogeneity of the medium is only in one direction, i.e. plane stratified 
plasmas (Stix and Swanson 1983, Swanson 1989, Stix 1992 , Budden 1985). The propagation of 
electromagnetic waves in four-dimensional space and time varying plasmas were studied before 
(Bernstein & Friedland 1983, Kravtsov 1969) in the context of geometrical optics. But in these 
studies the coupling of modes was not considered and the mere interest was concentrated an 
other themes like for example energy relations and the absorption of independent waves. There 
exist papers about four dimensional mode-conversion (Friedland, Goldner & Kaufman 1987, 
Kaufman & Friedland 1987), but the mathematical proof represented there is only valid for the 
very special cases where the coupling point is not a branch point. The reason is, that the different 
dispersion equations combined there result from one unique dispersion matrix, which is not the 
case when we consider for example the unified theory of a class of mode conversion problems, 
represented by Cairns & Lashmore-Davies 1983. Three dimensional transmission of the fast 
wave around the ion cyclotron resonance were considered by Friedland 1990 using reduction 
techniques. For recent investigations about combined mode conversion and absorbtion in one 
dimensional inhomogeneity see Cairns et.al. 1995 and references therein. 

If we want to study plasma physics more accurate, the propagation and coupling of electro- 
magnetic waves in media which vary in more than one dimension have many applications. The 
most important application of it is perhaps in radio-frequency heating and the investigations of 
plasmas in a tokamak (Cairns 1991). The plasmas there are considered as plane stratified, but 
being precise, we have also to take into account the curvature of the tokamak. In almost all 
cases, when wave coupling was under consideration, this curvature was assumed to be negligi- 
ble, or at most to give only small and unimportant effects. But experience in plasma physics 
has always shown, that small variations can give huge effects, because generally plasmas are 
unstable media. But also in the case of stability, it has been shown, that small variations can 
have essential new results, for example relativistic effects in wave absorptions (Mc Donalds et 
al 1994). Another point is, that in the papers about mode resonances, the magnetic field in 
the plasma is variable in only one direction. But we know, that in a okamak, the magnetic 
field is also toroidal and poloidal. Response functions for a variable toroidal magnetic field were 
calculated by Catto, Lashmore-Davies and Martin 1993. Another example of the application 
of wave coupling with general geometries is the ionosphere, which is usually considered as a 
plane stratified cold plasma (Budden 1985) and is a cost free medium for many plasma wave 
observations. Calculations taking into account the curvature of the atmosphere could give us 
new insight in the physics of the ionosphere. 
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It is the aim of this paper to give a systematic derivation of coupled wave equations in slowly 
varying plasmas which are three-dimensional inhomogeneous, non stationary, anisotropic and 
absorbing. In order to be self-contained, we include also a systematic representation of some 
results already presented on conferences before (Suchy & Sabzevari 1992 a,b) and refer to these 
papers for the explicit calculations. Some coupled wave equations were also represented there, 
but they must fulfil special conditions and are not applicable in general, as the derived equations 
in the present paper are. 



2 The Eikonal-Maxwell System of Differential Equations 

In a plasma with time and space dispersion the electric current density, dielectric displacement 
and magnetic induction must be written as integral equations 

J(X) = /VxV ker (X-X';X)-E(X') 

D(X) = /Vx' £ ker (X-X';X) -E(X') (1) 

B(X) = f d 4 X , n keT (X-X';X) -H(X') 

where the tensors of electric conductivity, dielectric - and magnetic permeability are written as 
the kernels of the integral equations. They have two kinds of dependencies on the space-time 
four-vectors X and X'. The ones before the semicolon are fast varying and are due to the 
collective effects of the medium, and those after the semicolon are slowly varying and due to the 
features of the medium itself, like density etc. In the integration over time, causality of course 
must be taken into account, but we redefine it in the kernels, so that the equations can be 
written in the above compact form, where all integrations are from — oo — > +oo. These equation 
must be solved together with the Maxwell equations 

VxH(X) = |D(X)+j(X) 
VxE(X) = -|B(X). 

To do this, we first transform the integral equation (|l]) into differential equations with the 
help of geometrical optics (Suchy & Sabzevari 1992a, Kravstov 1969), by writing j, D and B in 
eikonal form 

E(X) = E(;X)exp{#(X)} 

H(X) = H(;X)exp{i0(X)} (3) 
j(X) = j(;X)exp{i0(X)} 

where </>(X) is the eikonal, which in geometrical optics is a fast varying function of space and 
time, and E(; X), H(; X) and j(; X) are the amplitudes, which are slowly varying, because of the 
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slowness of variation of the medium. The first two equations of @ can also be considered as 
the wave entering the plasma, and the last as the movement of the charged particles in reaction 
of the waves. On principle, constant external fields could be added to (|3|), but later they give 
us no new physical results. 

Three different kinds of scales are considered here. First, we suppose that on each particle only 
the particles in the vicinity of that particle in space and time have an essential influence. This 
scale we define by Ai, which is small compared to the size of the plasma in space and time. 
Hence in (||) 

X-X' = 0(Ai)<l. (4) 

Then, it is sufficient to integrate only over the vicinity of the point X. The second scale is due 
to the slowness of the variation of the medium, which we denote by a small number A2. So we 
can write 

ggQX) gHQX) gGQX) 

-dx~ ' -ax~ ' ~ax- = 0(A2) « 1 (5) 

where G(;X) is an arbitrary function of the medium. The third scale, which we denote by 
A, is due to the largeness of the frequency and wave number, which are defined as the four 
dimensional wave vector 

K(X, = MS. (6, 
In geometrical optics, this can be written as 

K(X) = AG(;X) A>1 (7) 

where G(; X) is a slowly varying function due to the slowness of variation of the medium and is 
of order 0(1). From (0) and (|7|) 

0(X) = J K(X') • dXf = A J G(; X) • dX' » 0(A). (8) 
Hence, because of @ 

A K (X) = A^G(;X) = 0(AA 2 ). (9) 
In general, these three scales do not have to be of the same order 

0(Ai) + 0(A 2 ) ± 0(1). (10) 

Especially, the variation of the wave vector does not have to be as small along the ray path, as 
often is believed and assumed. 

Expanding now 0(X'),E(X') and H(X') around the point X, we obtain 

E(X) = e #(x) e i(X'-X)-K(x) | E , ; x) + / X / _ x) . 9 

I dX 

% - [(X - X)(X - X) : A K ] E(;X) + 0(Ax 2 A 2 ) + 0(AA 3 A 2 2 )| , 
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where we have used dyadic and double scalar products. A similar expression can be obtained 
for H(X'). Fourier transforming the fast varying part of <r^ er (X — X'; X), e^ er (X — X'; X) and 
//^ er (X— X'; X) in ([!]) and using the features of the delta function <5(K' — K) and it's derivatives, 
we can transform (Q) into (Suchy & Sabzevari 1992a, Kravtsov 1969) 

j(X) = e ! W[,(K;X)(l-ii-^)].E(;X) 

D(X) = eW[ £ (K;X)(l-ii-^)].E(;X) (12) 

B(X) = e W^(K;X)(l-ia~S)].H(;X) 

where 0(Ai 2 A2) + 0(AA 3 A2 2 )} terms have been neglected and the Janus- and Spreading- 
Operators are defined correspondingly as 

<-» d d <- d 2 d 

To combine them with the Maxwell equations, we first have to calculate the time derivatives 
of D(X) and B(X). To do this, we first write D(X) (and correspondingly B(X)) in the form 
(Suchy & Sabzevari 1992 a, Kravtsov 1969) 

D = D(K; X) exp{i<£(X)} (14) 

with 

D(K; X) = [e(K; X)(l - i d ~d)] ■ E(; X) + 0(Ax 2 A 2 ) + 0(AA! 3 A 2 2 ) (15) 
From (0), we can write 

m- ' 16 > 

" ' O(JL), (17, 



dKdK A 2 dGdG 'A 2 

and from the last two equations follows 

|D(K; X) = I [ e (K; X) • E(; X)] + O(^). (18) 

In all equations until now, we kept terms of order 0(Ai_A2), i.e. terms which are products of 
two small numbers. The last term in (]l~8|) is also the product of two small numbers. But we 
can neglect it compared with the last one, since 0(1/ A) -C O(Ai). This follows from the fact 
that in geometrical optics the scale of the wavelenght and period can be considered as much 
smaller than the scale of the distance and time, on which different parts of the medium have an 
essential influence on each other. Using the result ( |i~8| ) and corresponding expressions for the 
magnetic induction, we finally obtain 

J^D(X) = exp{i0(X)}e(K; X) { -iw(X)e(K; X) (19) 



w(X)e(K; X)(d +|<9) + [|e(K; X)] k } • E(; X 
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and a corresponding expression for <9B(X)/<9i, where e is interchanged with [i and E(;X) with 
H(;X). The time derivatives of D(X) and B(X) together with j(X) from (|l2|) can now be 
combined with the Maxwell equations (||) to give a system a partial differential equation for 
the components of the six dimensional electromagnetic field amplitude 



f(;X) 



E(;X) 
H(;X) 



(20) 



Since we have different modes propagating in the medium, we first expand the electromagnetic 
field into these modes 

f(X)=5>« X )f a (;X) (21) 

a 

where c/> a (X) is the eikonal and f a (;X) the electromagnetic amplitude for each mode and each 
mode has a separate wave vector 



K Q (X) 



«(X) 



(22) 



Then the system of partial differential equations for the electromagnetic amplitudes of the modes 
can be represented as 



1 



]T exp{# a (X)} M(K a ; X)(*+ d a +-d, 



0_ 

Of 



C(K a ;X) 



K 



f Q (;X)=0 (23) 



which we call the Eikonal-Maxwell system (Suchy & Sabzevari 1992a), with the Maxwell tensor 



M(K Q ;X) 



V 



u; a (k;X)e(K Q ;X)+i<T(K Q ;X) 



k x I 



k x I 



(24) 



w a (k;X)^(K a ;X) j 



where I is the unit tensor and k is the three dimensional spatial part of the wave vector (|2 
and w Q the frequency Q C(K a ; X) is the material tensor 



C(K a ;X) 



/ s(K a -X) \ 



^ ^(K a ;X) j 



(25) 



d a and d a are the operators ( |l3| ) corresponding to each mode. The derivatives to each mode 
are defined as 

Ft r) 

(26) 



^M(K a ;X) = AM(K;X) 



K=K Q 



1 In K Q , it is sufiicent to denote the index a to one component only, because of the dispersion relation (B 
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3 The eigenvalue spectrum of the modes and polarization vec- 
tors 



It can be shown (Suchy & Sbzevari 1992a) that the Maxwell- Tensor ( g4j ) can be written as the 
sum of two parts 

M(K;X) =M T (K T ;X) - AM L (K T ;X) (27) 
where we habe devided the wave- vector into a longitudinal and transversal part 

K = K T + Ag A (28) 

where g A is the axis in the direction of the longitudinal part. For example, A can be chosen as 
minus the frequency uj. Then Kj- is the spatial wave vector k and we call these frequency modes. 
For the case of stratified (not necessarily plane) media A can be chosen as the component of k 
perpendicular to the stratification surfaces. For a thorough discussion of various possibilities see 
Suchy & Sabzevari 1992a. It should be noted, that using the method represented here, restricts 
the dispersion of e, a and \i to the three components of JZt- 



The form ( p7| ) of the Maxwell tensor is called a generalized characteristic matrix. We can define 
right and left eigenvectors 

M(K a ;X)-g Q = [M T (K T ;X)-A Q (K T ;X)M L (K T ;X)]-g Q =0 

(29) 

g a -M(K a ;X) = ga -[M T (K T ;X)-A a (K T ;X)M L (K T ;X)]=0. 



Equations (|29| ) have nontrivial solutions only when 

DetM(K a ;X) = 0. (30) 

This gives us the dispersion relation. We choose the g a as the polarization vectors. The advan- 
tage of this choose is, that the right eigenvectors are parallel to the electromagnetic field in a 
homogeneous medium, as can be seen immediately from (23). We can then devide the space- 



time into small cells and in each cell the medium can be considered as homogeneous. When 
we go from one cell to the other, the g a will change direction and size. So the first equation 
of ( p9|) gives us the direction (but not the amount) of the wave field at each space-time point. 
The left eigenvectors have in general no real physical meaning. They are only useful for further 
calculations. But they can be considered as waves in a concomitant space (Suchy & Altman 
1975). 

It is easy to prove the following biorthogornality relation 

g a • M L • gp = 6 a p g a -M L - g a . (31) 
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The right and left eigenvectors are identical when the Maxwell tensor is equal to it's own trans- 
pose. In the case of a hermitean Maxwell tensor, the left eigenvectors become the complex 
conjugate of the right eigenvectors. 

Since the right (left) eigenvectors can be defined as an arbitrary linear combination of the 
columns (rows) of the adjoint of M(K Q ; X), we can define 



go 



AdjM(K a ;X)-c 
-v/c-AdjM(K a :X)-c 

(32) 

c-Adj M(K a ;X) 
-v/c-AdjM(K a ;X)-c 



where c and c are arbitrary constant column and row vectors. An expression for the biorthog- 
onality relation ( |3l| ) can be obtained (Suchy & Sabzevari 1992b) 

g Q -M L -g Q = 6 n(V-A 7 ). (33) 

The factor G = @(Ky;X) occures because we have a generalized characteristic matrix, i.e. the 
matrix (p7|). With this expression, it is possible to normalize the polarization vectors in a 
suitable form 



ga _ ga 



(34) 



and from the biorthogonality relations (31) follow biorthonormality relations 



f a • • f a = 1 always 

(35) 

f a -M L -f /3 = A Q / Xp. 
We define now the complete electromagnetic wave amplitude for each mode by 

i a = a a (;X) f a (K a ;X) (36) 

where a a is the scalar amplitude of the wave. To write the modes in the form of the last 
equation is of great advantage. The reason is, that f Q • ■ f a is equal to the energy of the wave 
propagating through the medium (Suchy & Sabzevari 1992a) when collisions are neglected, i.e. 
hermitean Maxwell tensor. Hence, because of the first of the biorthonormality rations in (|35|). 
|a Q | 2 is identical to the wave energy of the mode a. In the case of non-hermitean Maxwell tensor, 
this last quantity is proportional to the wave energy at each space-time point, hence |a a | 2 /|a#| 2 
gives us directly transmission-, reflection- or absortbtion- coefficients of different modes. 
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4 Mode coupling point 



Since the A's are functions of space and time, there are some points Xq, where A a (K^(Xo); Xo) = 
A / 3(Kr(Xo); Xo). Then from (p9|) and (^J) the polarization vectors f a and f« become parallel 
when the rank of the matrix M(K a ; X) (= M(Kg; X)) is equal to 5. This is because if the rank 
of a matrix is equal to the dimension minus the number of linearly independent eigencolumns 
of the matrix, then we have only one linearly independent f Q (= fg). Otherwise, there would 
be more than one linearly independent eigencolumns for each mode. For example, if the rank 
would be 4, then we would have two linearly independent eigencolumns for each a and /3, i. e. 
f«i (= f/3i) and f Q 2 (= f/32)- Hence f a i would be linearly independent of f^- In this case, at 
the point Xo, we can find two vectors f a and fg, corresponding to different modes, which are 
linearly independent, i.e. not parallel. No mode coupling would occure in this case. The two 
modes would pass by without any influence on each other (Budden 1985, chs. 17.6 and 17.7, 
Budden & Smith 1974). Hence the exact definition of a coupling point is 



Coupling point 44> {A a = \p for X = X and RankM(K Q ; X) = RankM(Kg; X) = 5} . 

(37) 



It should be noted, that at the coupling points the scalar amplitudes a a and ap are generally 
not equal. Physically, coupling does not occure at one mathematical point. Actually there is 
a small area around the coupling point, where coupling takes place. This coupling area can be 
defined as 



the coupling is weak. The more they get into the area, the stronger the coupling will be. At the 
coupling point, the coupling is maximum. Usually, a singularity is to be expected there (for the 
example of a stratified medium see Sabzevari 1992, 1993, Budden 1985, ch. 16.3, Budden 1972) 




From (|34|) and ( |35| ) follow the biorthonormality relations 



f«.M x -f/j = l for X = X 



(39) 



f a • M L ■ = for X/X . 



As we can see from (|39|), the point Xo is a "jump singularity", because the scalar product is 
discontinues there. 



2 work is in progress and will be represented in the framework of an arbitrary number of waves coupling to 
each other in a separate paper 
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5 Transport equations for the amplitude 



Since the polarization vectors f a can be calculated by the local eigenvector equation ( J29D at 
each space-time point X, we need another equation to calculate the scalar amplitudes a a along 



the rays to obtain the complete wavefield (36) at each space-time point. Putting (36) into the 
system (|23|), a system of partial differential equations for the amplitudes can be obtained (for 
the technical details of the calculations see Suchy & Sabzevari 1992a) 



. + r a 

P1 <9X p 7 



0. 



(40) 



The transport vectors W p7 give the change of the amplitude in the direction of this vector, as 
can be seen from the last equation itself. They have the form 



W 



(41) 



and the coefficients of the system are 



T n = 5 P1 ^^ + (X p - A 7 ) f p ■ M L ■ I f, 



d 



d 



dK T dK T ' <9X T 



k t +r 



in 



(42) 



with 



P7 



dr 7 y dr 7 y . 



M 



7 dK 



5 



+ 



/ \dx T "7 Kt 1 v dt j K 



dC\ 



(43) 



In these relations X = X^+gAX^ with gA-g = 1 (compare with (|28[)). r 7 is the parametrization 
of the ray 7 and is defined via the Hamiltonian equations (Suchy & Sabzevari 1992a, Berstein 
& Friedland 1983) 



where is defined as 



dr~„ 



dK 

dr~„ 



dD 7 (K;X) 
dK 



<9,D 7 (K;X) 



(44) 



ax 



K 



Det M(K; X) = -9(K T ; X) J| L> 7 (K; X) 



(45) 



with 



£> 7 (K;X) = A- A 7 (K T ;X). 



(46) 
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In the case where all the modes are propagating independently, i.e. when there is no mode 
coupling, there are no points where the ray paths coincide and the paths are independent from 
each other. In this case the functions 

7 (X) = / dX' • K 7 (X') (47) 

J ray path 

are independent from each other. The second terms in (f4l|) and ( |4"2| ) become identical to 
zero and the Kronecker symbol becomes one. From ([!(]) follows then an ordinary differential 
equation for the scalar amplitude of each mode a along it's ray path 

^ + T aa a a = (48) 



with 



The solution is 



8 D 

+ r aa . (49) 



a-a{T a ) = const exp < — / dT a T aa > (50) 



which gives us the amplitude along it's path. For more discussions about the solutions of 
independent modes see Suchy &; Sabzevari 1992a, sec. 7, Sabzevari 1992, 1993. 



Now the advantage of the representation ( pq ) (except that it is four-dimensional in space and 
time) becomes more transparent, since with the last equations, we can calculate the change 
of the scalar amplitude, and hence the energy directly, without bothering about polarization 
vectors 



6 Coupled wave equations 

In this section, we explicitly derive coupled ordinary differential equations for the amplitudes 
of coupled waves near their coupling point (except the coupling point itself, because there is a 
jump singularity, but we can approach the coupling point infmitly), i.e. where X a m \p. 



In the case of the coupling of the two modes a and /3, we can see from (47) that 4> a and (ftp can 



not be considered as independent anymore, because their ray path coincide in X- and K-space 
at the coupling point Xo- Since an arbitrary constant phase can be added to the eikonals, we 
can write near the coupling point in K-space 

0«(X) = / X dX' • K Q (X') « FdX! ■ Kp{X') = 4a(X). (51) 

J Xo J Xo 

As far as only the two modes couple and are decoupled of the other modes, <f) a and (frp are 
independent of the eikonals corresponding to the other modes. From (Eu) follows then 
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(52) 

W> a • §§f + • %g + r^ a a Q + T^a^ = 0. 
For the transport vectors follows from ( |4l"| ) 



W a/3 = (Ac - A/?) f a • M z • | t p -±- ) . (53) 

It may seem at the first moment, that this term is small, because A Q ~ \g. But this assumption 
is invalid, as we will see soon. To calculate the transport vector, we use (^) to write the 
equations as 

% = (A/3-A Q )~5 



fa — (A a — Xp) 

From (31) and ( |33| ) it is clear that 



(54) 

gg 



g a • M L • gp = for all A a and A/3 (55) 

We differentiat ip to Ky and multiply it from the right with f Q • Mi(Kr;X) to obtain with 
fl) for (H) 

W Q/3 = » — V 7 (56) 

©v/H-^,/? (^a ~ A 7 )(A/3 - A-y) 



Since gg is an arbitrary linear combination of the columns of AdjM(K / g; X), the elements of 
this column are polynomials of \a. Hence 

g ?dk = {w^) (HO + { g "dk) x • (57) 

It can be shown (see appendix) that near the coupling point in general 

dK T X a - \p 

This relation is not only valid for the derivatives to Ky, but for the derivatives of all kind of 
variables. The coefficients of the polynomials in g^ can in general be assumed as analytical, 
since they depend only on the elements of the matrices Mj(Ky;X) and Mi(Kr;X), which 
depend only on the e, fi, a, Ky and X. Therefore the first term in ([5?]) in general is growing 
relativ to the second as we approach the coupling point. Hence, we can write 

d ( d \ ( d\/3 



^"l^JteJ- (59) 
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Now from the first Hamiltonian equation (44) and the definition of Dp from ( |46"1 ) follows 

dX _ dXp 



drp 



+ gA 



(60) 



Near the coupling point (58) can be used to write (pO) as 



dX 



drp 8K.T 



(61) 



Putting this in (p9|), it follows from (|5 



W. 



a/3 ' 9X 



dap <y dap 

~7rt = 1 1 



where 



Ti 



To 



Vga-Mi-ga^-Mi-g^ 



-(A Q -A /3 )g )3 -M L .(^g Q ) 



w . _JL_ — jL 

Wq-q, ax — ^ 



(62) 



e \/ ri 7 ^ a , / 3( Aa_A ^)( A / 3 ~ A ^ 



(63) 



for A Q « Xp and A Q > (when X a ,< Xp, the sign of Ti and Y 2 would be interchanged). When 
only the two modes a and (5 couple, whilst the others are decoupled, Ti and Y 2 are finite. In 
contrast, by the same methods represented here, it can be shown that the T a p, Tp a , T aa and Tpp 
are growing to infinity as we approach the coupling point. This is due to the second derivative 



of the second term on the right hand side of (42). From (|5q ) this second derivative is of order 
0(A a — Xp)~ 2 and the singularity can not be removed when multiplied by (X a — Xp). 

We are now ready to construct the coupled equations. Near the coupling points the paths 
almost coincide and we can write dr a ~ drp = dr. Putting (]62] ) in (|52|), we obtain after some 
subtractions and additions the system 



d 

—a - T a 

or 







where 



with 
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r 



21 



I-T1T2 
1-TjTa 



12 



22 



1-T1T2 

1-TiTa ' 



(64) 



(65) 



(66) 
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The solution of (|64]) is 

a = exp|y rdr'j-ao (67) 
where ao is an arbitrary constant column. 

For a stratified medium, a different formalism than here, which was constructed for a coupling 
of an arbitrary number of coupled waves, also resulted in coupled wave equations with singular 
coefficients (Sabzevari 1992, 1993). But the form of the equation were of a kind, that the 
singularities first had to be removed in order to solve the problem (see also Budden 1985, 
ch. 16 and 17, Friedland 1985). This is not the case here in general. An example is given 
below. We believe, that the origin of these singularities, not only for the present case, but 
for all other cases studied before, could stem from the linearization of the coupling. Perhaps, 
these singularities would not appear, if the problem would be reconsidered under the aspect of 
nonlinearity (Tskhakaya 1996). These topics could be the basis for further investigations. The 
best method to investigate mode-conversion by the system of equations (^) is to use S-matrix 
methods. Using this method, there is no need to bother about the singularities at the coupling 
point. For an example see below. 



7 The example of a plane-stratified plasma 

To give a simple example of how to handle with the above formulas in practical cases and also 
to have a cheque of the correctness of the formalism, we consider the well known case of a plane 
stratified isotropic collisionless cold plasma with an arbitrary number of ion components. The 
magnetic field is in the z-direction normal to the stratification surfaces. The wave is in the 
x — z-plane and 9 the angle between the A;-vector and the z-direction. In this case, the Maxwell 
tensor has the simple form (Swanson 1989, ch. 2.1.2) 



oje P 











—k cos 9 








ue P 





k cos 9 





—A; sin 








uje P 





ksm9 








k cos 9 





UfiQ 








-k cos 9 





ksinO 





UflQ 








— ksin9 











U/JLq 



where 

3 

We consider monochromatic wave-guide modes, i.e. the component of the wave vector in the in- 
direction. Hence A = k z = k cos 9 and from Snell's law Kr = k x = ksin9 = (lj/co)sui#o =const., 
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where cq is the speed of light in vacuum and 9q the angle, the wave makes with the z-direction 
before entering the plasma, i.e. in the vacuum. Then, we have 



-1 





1 

-10 











(70) 



The coefficient (02n reduce to 



r P7 = f p .M L . (-f 7 



(71) 



To derive the right and left eigenvectors g and g, it suffices first to calculate the second row and 
column of Adj M respectivly. We define then 

" 
1 







c = (0,1,0,0,0,0) and c 



(72) 



and obtain from (32) 



go 



gc 



(P(P-nl)y- 



/^(P(P-nl))2n a cos I 



'^(P(P-n 2 a ))2n Q smt 



(73) 
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where n a = (co/uj)k a is the refractive index. From ([34]) we obtain 



(n a cos 9)~ 



co/^o 



lf Q (n aC os9) 



i _i 
j^na sin 9(cos 0) 2 



(74) 



Using (fr3), the equations (|63|) can be calculated 



+ n a rif3 



T 2 ± 2(n Q ^)V2 T P(P 



cos 2 # 



(75) 



and from (74) and ([71 



a,/3 



I I —{np cos 6). 



4^/n a np cos I np J dz 



(76) 



When we consider decoupled modes, i.e. cases where nowhere mode coupling occures, for the 
last equation follows F aa = Tpp = 0, and then from ( f48[) follows that the scalar part of the 
amplitude (but not the polarization vector) does not change when the wave propagates through 
the medium. This is to be expected, since wave absorbtion is a kinetic effect and does not occure 
in a cold plasma (see for example Antonson &; Manheimer 1978, sec. II), which also immediately 
can be concluded from the fact that the entropy can not change, when temperature is absent. 
In a cold plasma, mode coupling can occure only in the form of reflection and transmission, 
without absorbtion. 



We investigat now reflection at a cut-off. There, we have (Swanson 1989, ch. 2.1.2) 

n a = -na = P 1/2 -» 0. 
Hence (|64]) can be used, since (|5^) is valid. For (|75|) and (|76|) follows from (|77| ) 

2 = ±*f- 2 ' 



COS 
_ P" 



1 



f3a 



i d 
2PV2 ~dz 



ln(P 1/2 cos 



(77) 

(78) 
(79) 
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and from ©, (0) and © follows for (||) 

, / T, 1 \ / n„ \ 

= 0. (80) 




Using Snell's law, this can be written in the form of an equation for the difference of the energy 
of the incoming and reflected wave near the cut-off 

f] / \ _P 3 / 2 rip 

±ln(\a a \ 2 -\ap\ 2 ) = — T -?- for sinfl ± (81) 
dz V 2 sin 4 6*o 

From the definition of P in (p9l), it is clear that the derivative of P to z is always finite, hence 



-^ln(\a a \ 2 - | 1 2 ) -> for P -> 0. (82) 

(The case sin#o = can not be calculated by the above method, since in this case Snell's law 
is everywhere identical to zero. But physically, we can suppose 6q to be small and let the limit 
P — > go before 9q — > 0. Anyhow, this problem does not occure when calculating the S-matrix, 
which gives equivalent results, compare eq. fl85|)). 



As to be expected, the difference of the energy of the incoming and reflected wave is conserved 
at the cut-off, because there no part of the wave is transmitted. This is of course not the case, 
when we are away from the cut-off, where P is finite. There the difference of the energy of the 
incoming and reflected wave is not conserved anymore, because at each stratification surface a 
part of the wave is transmitted toward the cut-off. It is only at the cut-off, where pure reflection 
takes place. 

Solutions for the amplitudes a a and ap can be obtained near the coupling point by S-matrix 
methods (Volland 1962 a, b, Budden 1985, ch. 18.5). We define the cut off point by zq, i.e. 
P(zq) = 0. Z2 is a point above, and z\ a point below zq. The solutions near zq can then be 
written as 

a(z 2 ) = A(z 2 ,zi) -a(z{) (83) 
where the matrizant is defined as (Gantmacher 1960) 

rZi rZ2 fr\ 

A(z 2 , Zl ) = I+ drxT{Tx)+ dnVin)- dr 2 T(T 2 ) + ---. (84) 

./ 21 Jzx J Z\ 

Since zq is a cut-off, a a {z2) and ap{z 2 ) are the downgoing and upgoing modes above zq respec- 
tively and a a (zi) and ap(zi) the upgoing and downgoing modes below zq respectively (compare 
to Budden 1985, eq. 18.20). 
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For our purpose, it is sufficient to calculate the matrizant up to first order. We obtain 



A(z 2 ,zi] 



A(z 2 ,zi) 



P 5 / 2 {z 2 )-P 5 / 2 (z 1 ) -i{P 9 / 2 {z 2 )-P 9 / 2 (z 1 )) 



10 sin 4 



i(p9/2 (22) _p9/2 (2l)) 
18 sin B O 



18 sin 6 i 



-. _ P 5 / 2 (z 2 )~P 5 / 2 (z 1 ) 
10 sin 4 d 



pl/2( Z2 )_pl/2( Zl ) -j(p3/2( Z2 )-p3/2( 2l) ) 



j(p3/2( 22) _p3/2 (2i)) 

6 



1 



P 1 / 2 (z 2 )~P 1 / 2 (z 1 ) 
2 



for sin 9q / 



for sin 6q = 



(85) 



The diagonal elements of the matrizant are the reflection and the off-diagonal elements the 
transmission coefficients (Budden 1985, eq. 18.21, Volland 1962 a, eq. 21). We see, that at the 
cut-off, where z±,z 2 — > Zq, i.e. P — > 0, the reflection coefficients become 1 and the transmission 
coefficients 0, as is to be expected. 

Appendix 



Equation (45) can be written with the help of (|4"6|) as 



Det M(K; X) = 6 (A — A a )(A — \p) J[ (A - A 7 ) = Q(X 2 + BX + C) \ } (A - A n 

where B = S(Kr;X) and C = C(Kj , ;X) are analytical. The roots A a and A^ are 

-B ± v 7 ^ 2 - 4C 







At the coupling point, we have 



B 2 - AC = 



(86) 



(87) 



Hence the coupling point is a branch point (for the example of a plane stratified medium see 
Budden 1985, sec. 16.3 or Budden 1972, Budden & Smith 1974). We have 



8 1 dB 1 

dK T p ~ ~2 dK T ± A 



a 



8Kn 



(B 2 - AC) 



1 



A a — X/3 



(89) 



from which near the coupling point follows equation (|5 
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